COMPLEX ORGANIC MOLECULES AT HIGH SPATIAL 
RESOLUTION TOWARD ORION-KL II: KINEMATICS 
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ABSTRACT 



It has recently been suggested that chemical processing can shape the spa- 
tial distributions of complex molecules in the Orion-KL region and lead to the 
nitrogen-oxygen "chemical differentiation" seen in previous observations of this 
source. Orion-KL is a very dynamic region, and it is therefore also possible that 
physical conditions can shape the molecular distributions in this source. Only 
high spatial resolution observations can provide the information needed to disen- 
tangle these effects. Here we present millimeter imaging studies of Orion-KL at 
various beam sizes using the Combined Array for Research in Millimeter-Wave 
Astronomy (CARMA). We compare molecular images with high spatial resolution 
images that trace the temperature, continuum column density, and kinematics of 
the source in order to investigate the effects of physical conditions on molecular 
distributions. These observations were conducted at A = 3 mm and included 
transitions of ethyl cyanide [C2H5CN], methyl formate [HCOOCH3], formic acid 
[HCOOH], acetone [(CH 3 ) 2 CO], SiO, and methanol [CH 3 OH]. We find differences 
in the molecular distributions as a function of each of these factors. These results 
O . indicate that acetone may be produced by chemical processing and is robust to 

large changes in physical conditions, while formic acid is readily destroyed by 
gas-phase processing in warm and dense regions. We also find that while the 
spatial distributions of ethyl cyanide and methyl formate are not distinct as is 
suggested by the concept of " chemical differentiation" , local physical conditions 
shape the small-scale emission structure for these species. 

Subject headings: astrochemistry — ISM:individual(Orion-KL) — ISM:molecules — 
radio lines:ISM 
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1. Introduction 

The Orion-KL region is a very dynamic system with complicated molecular distributions, 
complex velocity structure, and extreme temperature and density gradients. Orion-KL has 
traditionally been viewed as a source that displays a high degree of nitrogen-oxygen "chem- 
ical differentiation" where oxygen-bearing molecules like methyl formate [HCOOCH 3 ] and 
nitrogen-bearing molecules like ethyl cyanide [C2H5CN] are spatially distinct. It has been 
suggested that this "c hemical different iation" can be driven by the molecular processes at 



work in these regions (jNeill et al.ll201ll ). with gas-phase reactions potentially accounting for 
the spatial separation of methyl formate and formic acid [HCOOH] emission peaks in this 
source. However, recen t high-resolution observations of molecular emission in this region 



f lFriedel fc Snvderl 120081 : iFriedel fc Widicus Weaver! 1201 ll ; IWidicus Weaver fc Friedell 120121 ) 



have called into question the traditional view of "chemical diff erentiation" in this source 



The m ost recent set of high spatial- resolution images presented by lWidicus Weaver fc Friedel 



( 120121 ) (hereafter Paper I) indicate that while the emission peaks for O-bearing and N-bearing 
molecules are spatially-separated in Orion-KL, the bulk of the complex molecular emission 
is co-spatial; and that any spatial differences are limited to smaller, localized regions. It is 
unclear from these observations whether chemical or physical processes shape these molec- 
ular distributions. While gas-phase chemical processing is a likely factor, it is also possible 
that varying physical conditions, including temperature, density, and kinematics, could in- 
fluence the molecular spatial distributions in these regions. However, only a limited amount 
of high spatial resolution information is available for this source, making definitive conclu- 
sions regarding the relative influence of these factors extremely difficult to determine. High 
resolution observations that probe the physical influence on the chemistry are necessary to 
disentangle and interpret these complicated processes. 

Here we extend the high spatial resolution observations of Orion-KL reported in Paper I 
to examine the physical conditions of this region. Using observations of methanol [CH3OH], 
we probe the temperature and velocity structure of this region. We also use the velocity 
information gained from these observations to develop a kinematic model for Orion-KL, and 
examine the influence of these kinematics on the molecular distributions in this source. We 
then compare this information with the molecular emission images reported in Paper I. From 
these results, we draw conclusions as to the mechanisms that are shaping the spatial distri- 
butions for each molecule observed. The observations, results, and discussion are presented 
below. 
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2. Observations 

The observations were conducted with the Combined Array for Research in Millimeter 
wave Astronomy (CARMA), and the details of the observations are presented in Paper I. In 
addition to those observations, we conducted observations in the CARMA C configuration 
in 2010 April. These observations were comprised of three 4.5 hour tracks. The observations 
have a phase center of a(J2000) = 05 /l 35 m 14 s .35 and 5(J2000) = -05°22'35".0 and have a 
synthesized beam of ~ 2.2 x 2.0". The u — v coverage of the observations gives projected 
baselines of 3.9-91.0 kA (13-304 m). The correlator was configured to have six 31 MHz 
windows and six 125 MHz windows, all with channel spacings ~450 kHz. The antenna- 
based gain calibration was done by self-calibrating on the SiO maser in Source I at 86.243 
GHz. The solution was then bootstrapped to the other windows. Phase offsets between each 
band and the SiO band were calculated and removed by using observations of 0423-013. 
Uranus was used to calibrate the absolute amplitudes of 0423-013 and these are accurate 
to within ~20%. The internal noise source was used to correct the passbands of each 31 
MHz window, and observations of 0423-013 were used to correct the passbands of the 125 
MHz windows. The data were calibrated, c ontinuum subtracted, and imaged using the 



MI RIAD software package (ISault et al.lll995l ). Note that the flux of the point source BN 
(see iFriedel fc Widicus Weaver! ( 1201 ll )) is similar across all array configurations, indicating 
that the accuracy of the absolute flux calibration between the different tracks is better than 
95%. 

In order to characterize the temperature profile of this region, 13 methanol transi- 
tions corresponding to 11 separate spectral lines were observed simultaneously. Methanol is 
widespread throughout the Orion-KL region, and therefore is an ideal tracer for the physical 
properties of the region. These 13 transitions span 2.3 orders of magnitude in upper state 
energy, thus providing a broad range of samples. Calibration uncertainties between the lines 
were minimized by conducting all methanol observations simultaneously, leaving only the 
true emission differences as differentiating factors. Tabled] lists the observed methanol tran- 
sitions, and includes the rest frequencies, quantum numbers, symmetry label, upper state 
energy, and linestrength information for each transition. 
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Table 1. Methanol Transition Parameters 



Quantum 


Frequency 


E u 


S^ 2 


Numbers 


(MHz) 


(K) 


(/J 2 ) 


19-2,18 " 192,17 


87,066.819 (10) 


467.7 


0.20 


18l0,8 - 18ll,7 


87,094.825 (5300) 


1308.0 


8.98 


13l0,3 " 13n,2 


84,930.9460 (10) 


1123.3 


4.01 


2-1,2 - 1-1,1 


96,739.363 (5) 


12.5 


1.21 


2o,2 - lo,i A+ 


96,741.377 (5) 


7.0 


1.62 


2o,2 - lo,i 


96,744.549 (5) 


20.1 


1.62 


2i,i - li,o 


96,755.507 (5) 


28.0 


1.24 


83,5 " 92,7 


94,541.806 (5) 


131.3 


2.24 


2i, 2 - li,i A+ 


95,914.310 (5) 


21.4 


1.21 


19 7 ,i3 - 20 6il4 A+ 


94,815.075 (50) 


684.8 


4.55 


197,12 - 206,15 


94,815.075 (50) 


684.8 


4.55 


246,19 - 237^6 


98,030.396 (118) 


888.5 


6.44 


24 6 ,i 8 - 23 7il7 A+ 


98,030.396 (118) 


888.5 


6.44 
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3. Results and Discussion 

The results of these observations were used to determine the effects of temperature, 
continuum column density, and kinematics on the molecular emission in the Orion-KL region. 
A temperature map of the region was constructed based on the methanol observations. The 
molecular emission maps presented in Paper I were then compared to this temperature 
map. The velocity information from the methanol observations was also used to examine 
the kinematics of the region. The results for each of these aspects of the observations are 
presented below. 



3.1. Rotational Temperatures 

Maps of CH 3 OH emission are shown in Figure [T] for the ^83,5-92,7 transition in i>lsr 
increments of 2.5 - 13 km s _1 . Assuming that the lines are optically thin and that the 
methanol is in local thermodynamic equilibrium, the total column density of CH3OH can be 
found from the relationships 
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T 
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where W is the integrated line intensity in Jy/beam km/s, Q rv is the rotational- vibrational 
partition function, E u is the upper state energy of the tra nsition in K, T r is the rotat ional 
temperature in K, C T is the o pacity correction fact or (see iGoldsmith fc Langerl (119991 ) ) . B 
is the beam filling factor (see lUlich fc Haasl (ll976l )). 6 a and 9 b are the FWHM Gaussian 
synthesized beam axes in arcseconds, Sfj, 2 is the product of the line strength and the square 
of the relevant dipole moment in D 2 , v is the transition frequency in GHz, and N u /g u is the 
upper state column density divided by its statistical weight (2J+1). 

A Boltzmann diagram was used to calculate the rotational temperature of the CH3OH 
transitions. Using the above equations, a plot of the natural log of Equation (JTJ versus E u 
was constructed. A weighted least squares fit to this plot yields the rotational temperature, 
which is equal to the negative inverse of the slop e of the line. The us e of this method does 
have some pitfalls which are described in detail by lSnyder et al.l (120051 ) . The most significant 
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Fig. 1. — J=8 3) 5-9 2 ,7 CH3OH maps for t> L sR of 2.5 - 13 km s 1 . Contours are ±3cr, ±6<r, 
±9cr,... (a = 194.5 mJy/beam). The synthesized beam is in the lower left corner of panel 
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factor to consider when using this analysis method is that the E u values cover a sufficient 
range to ensure a reliable temperature determination; this is not a concern in this work 
because of the large sample of lines observed. Also, the results of a Boltzmann analysis are 
more robust with a larger sample set. Therefore, in this work, the rotational temperature 
was only calculated in regions where at least 6 of the 11 lines were detected above the 3cr 
level. 

Figure H] shows the calculated rotational temperature for the Orion-KL region at rest 
velocities of 2.5 - 11.5 km s -1 . Figure |3] shows the uncertainty associated with these results. 
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Fig. 2. — Temperatures of the Orion-KL region measured from CH3OH observations. The 
bar on the right denotes the temperature scale in K. Regions which are white (black in online 
edition) did not have enough emission to make an accurate temperature determination. 



From these figures it can be seen that the CH 3 OH temperature measurements are not 
reliable toward the Hot Core, where the observational data were not of sufficient quality 
to enable a reliable temperature determination. This is because the lower E u transitions 
become highly self-absorbed in this dense region. Thus, while CH3OH appears to be a good 
tracer for temperature in most regions, it fails in the regions with highest methanol column 
density. In order to get an overall temperature view of the region, the ma ps from this 
work were combined with the NH3 temperature maps of iGoddi et al.l ( 2011b 1 ). The maps 
were combined in image space using the MIRIAD MATHS task ( Sault et al.l 1995 ). and the 
combined map is presented in Figure HI There is little spatial overlap between the two sets 
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Fig. 3. — Uncertainties of the calculated temperatures. 



of images except in the hot core region; where there was spatial overlap, the r esults agree to 



within the uncertainties. The peak temperature value determined was 415 K. iBeuther et al. 



(120051 ) also measured the temperature of the region with multiple CH 3 OH transitions at 
higher frequencies. This work more successfully determined the temperature in the hot core 
region, but suffered from opacity and double emission p eak issues in some of the more diffuse 
regions. The general temperature gradient observed by lBeuther et al.l (120051 ) and that found 
in this work are similar. The results presented here provide a more comprehensive overview 
of the temperature of Orion-KL. 
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Fig. 4. — T emperatures of the O rion-KL region measured from CH3OH and NH3 observations 
[NH3 from iGoddi et al.l (l2011bl )]. The bar on the right denotes the temperature scale in 
K. Regions which are white (dark blue in online edition) did not display enough emission 
intensity to make an accurate temperature determination. 
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3.1.1. Methyl Formate 

As with Paper I, all conclusions regarding HCOOCH3 apply equally to dimethyl ether 
[(CH 3 ) 2 0] , as they are co-spatial on all observed scales. Figure [5] shows the naturally weighted 
HCOOCH3 moment map (contours from Paper I) overlayed on the temperature map from 
above. From this map we can see that most of the HCOOCH3 emission arises from cooler, 
less dense regions such as IRc5 and IRc6. There is some emission arising from the warmest 
regions (Hot Core and IRc7); however, the peaks of this emission are not coincident with 
the temperature peaks. There is no temperature information corresponding to the most 
extended parts of the HCOOCH 3 emission near the Compact Ridge. This is likely due to 
these regions having low density, and the emission being extended and therefore resolved out 
in the CH3OH observations. 

3.1.2. Ethyl Cyanide 

As with Paper I, all conclusions regarding C2H5CN apply equally to vinyl cyanide 
[C2H3CN] and methyl cyanide [CH 3 CN], as they are all co-spatial on all observed scales. 
Figure E] shows the naturally weighted C2H5CN moment map (contours from Paper I) over- 
layed on the temperature map. Unlike HCOOCH3, C2H5CN peaks in the warmest, densest 
regions (Hot Core, SMA 1, and IRc7). The C 2 H 5 CN peak at ~12 km s" 1 (see Figure 10 of 
Paper I) is associated with the cooler region near IRc6. 

3.1.3. Acetone 

Figure [7] shows the naturally weighted (CH^CO moment map (contours from Paper 
I) overlayed on the temperature map. (CHa^CO appears to have some characteristics in 
common with both HCOOCH3 and C2H5CN, as it has strong peaks near the warmest regions, 
but also has some extended structure in the cooler regions. 

3.I.4. Formic Acid 

Figure |8] shows the HCOOH moment map (contours, from Paper I) overlayed on the 
temperature map. HCOOH has a unique distribution with regard to the temperature map 
compared to the other complex molecules observed. It has one weak peak near the warmer 
regions, but the bulk of the emission is coming from the cooler, more extended regions. Full 



05 h 35 m 15 s 14 s 8 1 4 S 6 1 4 S 4 1 4 S 2 1 4 S 1 3 S 8 1 3 S 6 

a (2000) 

Fig. 5. — Methyl formate [HCOOCH 3 ] contours overlayed on the color scale temperature 
profile of Orion-KL. The contours are 3a, 6a, 9a, ... (a=4.1 mJy/beam km s _1 ). The 
moment spans the entirety of the HCOOCH3 emission in velocity space (1.9 - 11.8 km s _1 ). 
The color bar on the right is the key for the temperature scale in K. 
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Fig. 6. — Ethyl cyanide [C2H5CN] contours overlayed on the color scale temperature profile 
of Orion-KL. The contours are 3a, 6a, 9a, ... (a=4.1 mJy/beam km s _1 ). The moment 
spans the entirety of the C2H5CN emission in velocity space (-12.8 - 18.4 km s _1 ). The color 
bar on the right is the key for the temperature scale in K. 



13 




05 h 35 m 15 s 14 s 8 1 4 S 6 1 4 S 4 1 4 S 2 1 4 S 1 3 S 8 1 3 S 6 

a (2000) 



Fig. 7. — Acetone [(CH 3 ) 2 CO] contours overlayed on the color scale temperature profile of 
Orion-KL. The contours are 3a, 6a, 9a, ... (a— 8.1 mJy/beam km s _1 ). The moment spans 
the entirety of the (CH^CO emission in velocity space (2.4 - 8.7 km s _1 ). The color bar on 
the right is the key for the temperature scale in K. 
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temperature measurements are not available for much of the HCOOH distribution because 
large structures were resolved out in the CH 3 OH observations. 



3.1.5. SiO 



Figure M shows the moment map SiO contours overlayed on the temperature map. The 
SiO distribution peaks adjacent to, but not co-spatial w ith, the higher temperature and 
continuum column density regions. iGoddi et all ( 1201 lbl ) conclude that the main heating 
source in the Orion-KL region is shock heating. Thus it is quite likely that the SiO is tracing 
the shocks that heat the region as they interact with the dense ambient material. 
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Fig. 8. — Formic acid [HCOOH] contours overlayed on the color scale temperature profile of 
Orion-KL. The contours are 3a, 6a, 9a, ... (<r=8.1 mJy/beam km s _1 ). The moment spans 
the entirety of the HCOOH emission in velocity space (0.2 - 8.1 km s _1 ). The color bar on 
the right is the key for the temperature scale in K. 
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Fig. 9. — SiO contours overlayed on the color scale temperature profile of Orion-KL. The 
contours are logarithmically scaled in order to show all of the weaker details. The moment 
spans the entirety of the SiO emission in velocity space (-43.9 - 54.9 km s _1 ). The color bar 
on the right is the key for the temperature scale in K. 
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3.2. Kinematics 

These new observations enable discussion of the kinematics of the Orion-KL region. As 
SiO is a well-known shock tracer, we compare the distribution of each molecule with that of 
SiO to gain insight into the physical conditions for the regions where the different molecules 
are detected. 



3.2.1. Methyl Formate 

Figure [10] shows the naturally weighted moment map of HCOOCH3 emission overlayed 
on the moment map of SiO emission. The HCOOCH3 emission appears to be coming from 
the edges of the dense SiO outflow. In these regions, the outflow is impacting the ambient 
dust and gas, and is likely releasing the HCOOCH3 molecules from icy grains and into the 
gas phase. 



3.2.2. Ethyl Cyanide 



Figure [TTJ shows the C2H5CN moment map overlayed with the SiO outflow moment 
map. From this map, it appears that the C2H5CN emission is coming from the regions 
where the outflow is impacting the ambient dusty material to the northeast of Source I and 
near IRc7, releasing the C2H5CN into the gas phase. The lack of detection of C2H5CN to 
the west /northwest of the outflow (mirroring what is seen to the east /northeast) can be 
explained by a lack of dusty material in this direction (see Section [3.3. 51) . This hypothesis 
is supported by the fact that I Rc3 appears to be a r eflection of IRc2, indicating little-to-no 
dust in the intervening regions (jSimpson et all 120061 ). 



3.2.3. Acetone 

Figure [121 shows the moment map of (CH 3 ) 2 CO emission overlayed on the SiO outflow 
moment map. Although the (CH 3 ) 2 CO image has a much lower spatial resolution than 
C2H5CN and HCOOCH3, it is still clear that it is coming from the regions in the center of 
the outflow, corresponding to the densest regions where nearly all of the molecular emission 
features observed in these studies overlap. However, (CH 3 ) 2 CO is not present in all of the 
dense regions, while C 2 H 5 CN is present in all of these regions. (CH 3 ) 2 CO does not appear 
to have any significant emission to the northeast of Source I. 
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a (2000) 

Fig. 10. — Naturally weighted, methyl formate [HCOOCH3] emission (black contours) over- 
layed on SiO emission (grey contours, red contours in online edition). 

3.2.4- Formic Acid 

Figure H3] shows the HCOOH moment map overlayed on the SiO outflow moment map. 
While the HCOOH emission appears near the edges of the outflow, its emission is notably 




Fig. 11. — Naturally weighted, ethyl cyanide [C2H5CN] emission (black contours) overlayed 
on SiO emission (grey contours, red contours in online edition). 



farther away from the center of the outflow than the emission regions corresponding to the 
other molecules. These regions traced by the HCOOH emission may still contain material 
from the outflow, but at a low enough density as to be undetected in these observations. 
This indicates that HCOOH may be easily liberated from grain surfaces by shocks, but is 




Fig. 12. — Naturally weighted, acetone [(CH^CO] emission (black contours) overlayed on 
SiO emission (grey contours, red contours in online edition). 



then quickly processed into other molecules, as indicated by the lack of HCOOH emission 
from behind the initia l shock front. This is chemically consistent with the results of recent 
astrochemical models ( Garrod. Widicus Weaver. Sz Herbst II2OO8 : Laas et al. 2011 ) and with 
recent observational studies (INeill et al.ll201ll ). which suggest that HCOOH is highly reactive 




Fig. 13. — Formic acid [HCOOH] emission (black contours) overlayed on SiO emission (grey 
contours, red contours in online edition). 
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3.3. Continuum 



Comparison of molecular distributions with continuum sources can also give insight into 
the physical conditions of the different regions. In this section we compare the molecular 
distributions o f the observed species with t hat of the A = 3 mm continuum. The continuum 



maps are from lFriedel &: Widicus Weaver! (120111 ) 



3.3.1. Methyl Formate 

Figure [H] shows the naturally weighted HCOOCH3 contours overlayed on the gray scale 
A = 3 mm continuum. This map shows that the HCOOCH3 emission peaks are closely 
associated with the weaker continuum peaks, but not the strongest ones. 



3.3.2. Ethyl Cyanide 

Figure [15] shows the naturally weighted C2H5CN contours overlayed on the gray scale 
A = 3 mm continuum. This maps shows that the C2H5CN peaks are nearby to strong 
continuum sources. 



3.3.3. Acetone 

Figure [TBI shows the naturally weighted (CH 3 ) 2 CO contours overlayed on the gray scale 
A = 3 mm continuum. This map shows that (CH 3 ) 2 CO peaks near the continuum, but not 
with the continuum peak. 



3.3.4. Formic Acid 

Figure [T71 shows the HCOOH contours overlayed on the grayscale A = 3 mm continuum. 
The map shows that there are a few weak HCOOH peaks associated with the more extended 
continuum, southwest of Source I and near IRc3. However the bulk of the HCOOH emission 
is not associated with any detected continuum emission. This is not an issue of resolution, 
as some of the continuum data were obtained at the same time as the HCOOH data. 
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Fig. 15. — Naturally weighted C 2 H 5 CN contours overlayed on the naturally weighted A = 3 
mm continuum gray scale. 
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Fig. 16. — Naturally weighted (CH 3 ) 2 CO contours overlayed on the naturally weighted A = 
3 mm continuum gray scale. 



-26 - 




Fig. 17. — Naturally weighted HCOOH contours overlayed on the naturally weighted A = 3 
mm continuum gray scale. 



-27- 



3.3.5. Global Kinematics 



It is apparent from previous observations that the Source I and BN sources within the 
Orion-KL region had some type of dyn amic interaction ^500 years ago and are moving 
apart at ~0.02" (~50 km s~ x ) per year (IPlambeck et al.lll995l ; iGoddi et al.l l2011al ) . There 
are also two known outflows in the Orion-KL complex, both centered near Source I. The first, 
the "high velocity" (30-200 km s _1 ) outflow, is traced by CO, and possi bly SiO, and has a 
very l arge opening angle (~1 rad) in the northw est /southeast directions ( iChernin fc Wright 



19961 ; IZapata et al.l l2009t IPlambeck et al.l 120091 ). The second, th e "low velocity" outfl ow, is 
perpendicular to the high velocity outflow and is traced by SiO (iGreenhill et al.lll998l ) . The 
actual velocity of the low velocity outflow is not well-co nstrained. IGreenhill et al.l (120041 ) 



concluded the outflow velocity was ~18 km s 1 , however Plambeck et al. ( )2009l ) concluded 



that the outflow velocity must be much larger else there would be a "trail" from the SiO 
outflow tracing the proper motion of Source I. IPlambeck et al.l ( 120091 ) also suggested that 
the two outflows may be related as the SiO outflow may be precessing rather quickly, and 
the high velocity outflow may be from an era when the SiO outflow was oriented nearly 
perpendicular to its current position. We hypothesize that there is yet another interpretation 
of these high velocity clumps. These clumps might have had close encounters with regions 
of higher density/mass that may have altered their trajectory, both within the plane of the 
sky and in/out of the plane of the sky. 

Figures [18] and [19] show a potential view of the general kinematics of the Orion-KL 



region . The gray scale is a moment map of the SiO outflow [from Figure 3 of IPlambeck et al. 



( 120091 )]. and the black contours are a combined moment map of the C2H5CN and HCOOCH3 
naturally weighted data (Figure IT8|) and the A = 3 mm continuum (Figure IT9|) . The letters 
(A - E) indicate strong molecular peaks while the letters (a - e) indicate strong continuum 
peaks (excluding Source I and BN). The blue arrow and dashed lines i ndicate the most 
probable path of Source I after its interaction with BN (see IGoddi et al.l (l2011ah ) . and the 
blue star indicates the most probable location of this interaction. The solid red lines indicate 
potential paths taken by the SiO outflow. The broken red line indicates a potential path for 
a high velocity clump of gas after leaving Source I. Labels 1-4 indicate high velocity clumps 
(HVC). The outflowing gas in the HVC may encounter temperature/density enhancements 
as it interacts with the denser ambient clumps. This may be evident in our temperature map 
of the region (see Figure H|). We see no evidence of an enhancement for HVC 1-3, however 
the temperature maps do not cover most of the path covered by these HVC. For HVC 4, we 
do see a temperature enhancement leading from the hot core region southwest to the edge 
of the temperature map. The "bend" in this enhancement may be due to the transverse 
motion of Source I over time and may be tracing the history of travel for HVC 4. 



-28 - 



In Figure [TBI the northeast lobe of the SiO outflow seems to be impeded by the material 
near clump A, while much of the material seems to be "escaping" to the east of clump A. 
Additionally, HVC1 appears to have broken through the region between clumps A & B. 
Clump C appears to be an impediment to the outflow in the southwest while HVC4 seems 
to have passed the clump to the west. From Figure [191 we can see that the mass/density of 
clump 7 may be impeding the outflow while HVC4 passes to the west. 



3.4. Putting it All Together 

From these results, it is clear that the dynamics and chemistry of the Orion-KL region 
are very complex. Comparison of the molecular emission features with the kinematic tracers 
allows conclusions to be drawn about the effects of physical conditions on the chemistry of 
this region. Here we overview the major conclusions for each molecule studied. 



3.4-1- Ethyl Cyanide 

C2H5CN is present at highest abundance in the regions of highest temperature, along the 
denser edges of shocks, and near the higher column density dusty regions. The shocks are the 
most likely mechanism for releasing the C2H5CN from grain surfaces, as they provide more 
thermal and kinetic energy than any other source in the region. These results also indicate 
that C2H5CN must also forming through gas-phase processing of grain mantle material in 
the warmer regions of the source. Otherwise, C2H5CN would be destroyed in these warmer 
regions after being liberated from icy grains by shocks. 

In contrast is the emission from the high velocity component near IRc3/6/20. Here 
the temperature is lower, and this region is on the outermost shock region in an area of 
lower continuum column density. There is also little continuum emission here. These results 
suggest that there are several conditions under which C2H5CN can form and be released into 



the g as phase. This area may be indirectly heated by IR radiation from IRc2 (ISimpson et al. 



20061 ). Thus, it seems that C 2 H 5 CN can be present at a range of physical conditions, including 
very dense and hot regions and less dense and cooler regions. Regardless, a heating source 
seems to be required for its formation or liberation. 
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Fig. 18. — The gray scale is a moment map of the SiO outflow (from Figure 3 of 
Plambeck et al.l ( 120091 )) and the black contours are a combined moment map of the C2H5CN 
and HCOOCH3 naturally weighted data. The letters (A - E) indicate the main molecular 
peaks. The thin black (blue in online edition) arrow and dashed lines i ndicate the most 
probable path of Source I after its interaction with BN (see iGoddi et al.l (j2011al )). and the 
star indicates the most probable location of this interaction. The solid black (red in online 
edition) lines indicate potential paths taken by the SiO outflow, and labels 1-4 indicated 
high velocity clumps for reference in the text. The broken black (red in online edition) line 
indicates a potential path for a high velocity clump of gas after leaving Source I. 
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Fig. 19. — Th e gray scale is a moment map of the SiO outflow (from Figure 3 of 
Plambeck et al.l ( 120091 )) a nd the black contours are the n aturally weighted map of the A 
= 3 mm continuum from iFriedel fc Widicus Weaver! (1201 ll ). The Greek letters (a - e) in- 
dicate the stronger continuum peaks (excluding Source I and BN). The thin black (blue in 
online edition) arrow and dashed lines i ndicate the most probable path of Source I after its 
interaction with BN (see iGoddi et al.l (j2011al )). and the star indicates the most probable 
location of this interaction. The solid black (red in online edition) lines indicate paths were 
the SiO outflow appears to have broken through the surrounding material, and labels 1-4 
indicate high velocity clumps for reference in the text. The broken black (red in online 
edition) line indicates a potential path for a high velocity clump of gas after leaving Source 
I. 
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3.4-2. Methyl Formate 

The HCOOCH3 emission appears to be coming from two distinct regions. The first is 
quite extended, near the location of the Compact Ridge. This region is not associated with 
any strong continuum emission. This region is on the extreme leading edge of the outflow 
from Source I, indicating that the HCOOCH3 is being liberated from grain surfaces. This 
emission region is likely the primary contributor of flux in single dish observations, as it 
would be least affected by beam dilution. 

By comparing the compact HCOOCH3 emission with the different tracers, it appears 
that the vast majority of this emission comes from regions on the outer, lower column 
density portions of the SiO outflow; lower temperature (<150 K) areas; and areas of lower 
dust concentration. There is some notable emission from the hottest /densest parts of Orion- 
KL, but these regions coincide with the edge of the HCOOCH3 emission. It appears that 
HCOOCH3 is released into the gas phase by the initial expanding shock from the outflow, 
and is destroyed by the arrival of the higher density, higher temperature components of the 
outflow. HCOOCH3 is likely processed through additional gas-phase reactions once it is 
liberated from the grain mantle material. Unlike the case of C2H5CN emission peaking in 
the warmer regions, there are no signs of HCOOCH3 being formed in the warmer regions of 
the source, as its emission peaks are offset from these warmer regions and instead coincide 
with the leading edges of the shocks. 



3.4.3. Acetone 

While there are no corresponding high-resolution images for (CHs^CO, the existing 
images reveal that it does not have a distribution that matches either C2H5CN or HCOOCH3. 
Instead, acetone appears to peak in the regio n of most significant over lap between HCOOCH 3 



and C2H5CN (as was originally reported by iFriedel fc Snyder! (120081 )). The primary peak is 



located in the warm transition region between the hot and dense part of the outflow, and 
the cooler outer edges. However, there is notable emission from both the hottest and densest 
parts of the cloud, as well as the cooler and less dense regions. This indicates that (CH 3 ) 2 CO 
is present in a wide variety of physical conditions, and its emission may therefore be shaped 
primarily by chemical processes rather than the physical conditions of the cloud. 
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3.4-4- Formic Acid 



HCOOH has a unique distribution relative to the physical tracers when compared to 
the other molecules in this study. Its strongest emission peak is located on the extreme outer 
edge of the SiO outflow, where there is no detectable continuum emission. There are a few, 
weaker HCOOH peaks in the warmer denser regions, but overall it appears that HCOOH 
is easily released by the shocks from the outflow, and then easily destroyed in the higher 
density and warmer regions. This is consistent with che mical models that predict that formic 



acid i s readily processe d through gas-phase reactions (IGarrod. Widicus Weaver, fc Herbst 



20081 : lLaas et al. 1 120111 ). 



4. Conclusions 



We have presented CARMA observations of several molecular species, including ethyl 
cyanide, methyl formate, acetone, and formic acid, toward Orion-KL. These observations 
cover over an order of magnitude in spatial resolution, yielding excellent sensitivity to both 
small- and large-scale structures. We have compared these results with tracers of physical 
conditions (temperature, continuum column density, and dynamics). Overall, the results 
presented here indicate that the distributions of complex molecules trace the dynamics of 
the region, and that in some cases varying physical conditions can affect the morphology 
of molecular emission. These results indicate that all of the observed species can exist in a 
range of temperature and density conditions; HCOOH appears to have the narrowest range of 
stability in terms of temperature and density, while C2H5CN can exist in the widest range of 
conditions. Furthermore, it appears that chemical processing may produce (CH 3 ) 2 CO under 
a variety of physical conditions, while HCOOH is readily destroyed by gas-phase processing 
in warmer, dense regions. Lastly, it appears that the traditional view of nitrogen-oxygen 
"chemical differentiation" in Orion-KL, where different chemistries lead to spatially-distinct 
emission regions for these two classes of molecules, does not give the complete picture of this 
region. HCOOCH3 and C2H5CN are found to be primarily co-spatial in this source. Instead 
of being shaped by chemical processing, local physical conditions like temperature and den- 
sity most dramatically shape the relative distributions of these particular molecules. These 
results reveal that astrochemical modeling of sources must include the effects of physical 
conditions and differing spatial distributions before the observations of a given source can 
be compared directly to the model. Sophisticated hydrodynamic models will be required to 
fully interpret the results presented here. 
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